Dynamic generation of maximally entangled photon multiplets by adiabatic passage 



O 

o 
o 

Oh- 
< 

(N 



> 
00 

o 
^- 

O 
O 
O 

^ : 

Oh: 

^— > ' 
G ■ 
cd 

3 : 
cr. 



W. Langc* and H. J. Kimble 
Norman Bridge Laboratory of Physics 12-33, California Institute of Technology, Pasadena, California 91125 

(10 August, 1999) 

The adiabatic passage scheme for quantum state synthesis, in which atomic Zeeman coherences 
are mapped to photon states in an optical cavity, is extended to the general case of two degenerate 
cavity modes with orthogonal polarization. Analytical calculations of the dressed-state structure 
and Monte Carlo wave-function simulations of the system dynamics show that, for a suitably chosen 
cavity detuning, it is possible to generate states of photon multiplets that are maximally entangled 
in polarization. These states display nonclassical correlations of the type described by Greenberger, 
Home, and Zeilinger (GHZ). An experimental scheme to realize a GHZ measurement using coin- 
cidence detection of the photons escaping from the cavity is proposed. The correlations are found 
to originate in the dynamics of the adiabatic passage and persist even if cavity decay and GHZ 
state synthesis compete on the same time scale. Beyond entangled field states, it is also possible 
to generate entanglement between photons and the atom by using a different atomic transition and 
initial Zeeman state. 

PACS number(s): 42.50.Dv, 03.65.Bz 



I. INTRODUCTION 

A quantum system consisting of several components 
is called entangled if its state cannot be expressed as a 
direct product of substates for the constituents. Such a 
system has no classical analog and hence shows distinct 
quantum mechanical features. Of particular interest are 
multiplets of particles in maximally entangled states, in 
which a measurement of the energy eigenvalue of one par- 
ticle completely determines the outcome of the same mea- 
surement for all the others. An example in the case of two 
particles are the EPR or Bell states [1], which give the 
largest contradiction with the local realism assumption of 
Einstein, Podolsky, and Rosen [2]. Greenberger, Home, 
and Zeilinger [3,4] and later Mermin [5] have shown that 
a much stronger and more direct test of local realism may 
be achieved by using three or more maximally entangled 
particles. Violation of local realism can then be obtained 
in a single set of experiments, whereas EPR correlations 
between two particles show up only statistically in a se- 
ries of measurements. 

A number of different experimental schemes have been 
proposed to generate a maximally entangled or GHZ 
state involving three or more photons or atoms [3], Pro- 
posals for optical experiments typically use photon pairs 
generated in parametric down conversion to obtain GHZ 
correlations [6-8] or polarization correlations among pho- 
tons emitted in cascaded atomic decays [4]. Other 
schemes suggest utilizing the QED interaction of atoms 
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and field modes of a cavity to generate entanglement. 
Three or four two-level Rydberg atoms may be prepared 
in an entangled state through their interaction with the 
field in a microwave cavity [9-12]. Conversely, three mi- 
crowave cavities may be entangled due to the interaction 
with a Rydberg atom [13]. Yet another suggestion uses 
a single photon in a double Mach-Zehnder interferome- 
ter to create a three-atom GHZ state [14]. There are 
also mixed schemes with entanglement between different 
types of subsystems, for example four cavity modes and a 
single atom [15]. While in most proposals the entangled 
particles originate in a single localized source, GHZ cor- 
relations may also be produced using particles from sep- 
arate sources [16,17]. In this way recently photon states 
were realized that retrodictively reproduce GHZ statis- 
tics of photocounts [18]. Another system that should be 
capable of producing entanglement are two-level ions in a 
linear quadrupole trap, subjected to appropriately timed 
pulses [19]. Some of these schemes allow the subsystems 
to be probed in spatially separated locations, which is 
important if the issues of nonseparability and nonlocal- 
ity in quantum mechanics are to be investigated. In other 
cases such a separation of the subsystems is impossible 
in principle, thus precluding tests of the causal depen- 
dence between the measurements. In spite of the many 
theoretical proposals, however, to date no experimental 
realization of a GHZ state with three or more particles 
has been achieved. 

Here, we present a novel approach to the generation 
of GHZ correlations, adapting a method from the field 
of quantum-state synthesis in cavity QED [20-22]. The 
particles we consider are photons in a single longitudinal 
mode of an optical resonator. One way to create "arbi- 
trary" nonclassical states of the field in such a system 
is the method of adiabatic passage [20,23,24]. In this 
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protocol, an atom that is prepared in a coherent super- 
position of Zeeman ground states sequentially traverses 
a resonant quantized cavity mode and a coherent-state 
laser field with overlapping spatial profiles. If the transit 
is slow enough (adiabatic), the system remains in a sin- 
gle quantum state and the coherence of the Zeeman sub- 
levels is mapped to the photon distribution of the cavity 
mode. Therefore, the number of quantum states of the 
field that may be synthesized is only limited by the pos- 
sible Zeeman state superpositions of the atom. Previous 
work [20,23] shows this scheme to be robust, without the 
requirement for precise control of pulse area or timing. 

In our scheme for the generation of a GHZ state, the 
cavity photons are entangled in polarization. Since in the 
previous theory of adiabatic passage [20,23] only a sin- 
gle polarization was considered (which was taken to be 
cr_, implying an infinitely detuned <7 + -polarized mode), 
we start in Sec. II by extending the method to the case 
of two degenerate orthogonal polarizations for the cavity 
mode. The model allows us to find a straightforward pro- 
cedure for the synthesis of a GHZ state of polarization 
entangled photons and to calculate the efficiency of its 
generation. Throughout the paper we utilize parameters 
already realized in present day cavity QED experiments. 
Note that from an experimental point of view, our re- 
sults are relevant even for the generation of Fock states 
of a single polarization, since in realistic Fabry-Perot res- 
onators, modes with both transverse polarizations are al- 
ways present. Although in practice small biases towards 
one polarization may occur (e.g., due to birefringence of 
the mirrors), our model assumes that the frequency split- 
ting of the polarization modes of the cavity is negligible 
(degenerate case). However, in Sec. IV we present numer- 
ical calculations for finite mode splitting (in the circular 
polarization basis) to derive a quantitative criterion for 
the validity of the degenerate model. These results also 
show how large a mode splitting is necessary for a single 
transverse mode description of the system to be applica- 
ble. 

In the single polarization case [20,23], the adiabatic dy- 
namics evolves in a one-dimensional space along the so- 
called dark state, a superposition of ground states which 
is decoupled from excited atomic levels. Each atomic 
ground state within the manifold of Zeeman substates 
is correlated with one (and only one) photon number of 
the cavity field, as any given Zeeman state is reached 
by a unique sequence of photon absorption and emission 
cycles. Thus a shift of the atomic population is always 
accompanied by a corresponding shift of the photon dis- 
tribution. By contrast, with both (degenerate) polar- 
izations of the cavity field taken into account, the sys- 
tem behavior is considerably more complex. Population 
may be shifted to higher or lower magnetic sublevels, de- 
pending on the polarization of the photons involved. As 
a consequence there is no intrinsic limit to the number 
of photons that may be generated by adiabatic passage. 



The resulting final states of the cavity field are no longer 
expected to be pure Fock states, but a complex superpo- 
sition of number states for both polarizations, even when 
the initial atomic state is a single Zeeman substate. This 
is shown in Sec. Ill by analyzing the structure of the 
relevant eigenstates. 

After characterizing the dynamics of state synthesis by 
adiabatic passage with two polarizations, we proceed by 
proposing a new strategy to arrive at pure Fock states 
of the cavity field by making use of a suitable detun- 
ing of the cavity field. The feasibility of the method is 
demonstrated in Sec. IV by a quantum Monte Carlo sim- 
ulation of the system dynamics. An important result is 
that there is no restriction on our method with respect to 
polarization: the scheme is capable of generating num- 
ber states of either cavity polarization or even arbitrary 
superpositions of both. The last fact opens the door for 
the generation of a GHZ state of polarization entangled 
photons starting from a suitable superposition of Zeeman 
ground states. Details of this new protocol are explained 
in Sec. V. 

In a realistic experiment, the synthesis of a GHZ state 
in a cavity will certainly be affected by energy dissipation. 
In Sec. VI we incorporate this fact into the model and 
utilize it to detect GHZ correlations by coincidence mea- 
surements of the photons escaping from a leaky cavity. 
Somewhat surprisingly, we find that even when decay sets 
in before the generation of the GHZ state is finished (i.e., 
when the cavity decay and adiabatic passage times are 
comparable), correlations between the photons persist. 
Obviously, these correlations originate in the dynamics 
of the adiabatic process rather than in the properties of 
a stationary final state of the system. 

A variation of the method that creates entanglement 
among photons and atoms is presented in Sec. VII. Here, 
the initial state for the adiabatic passage is a single Zee- 
man ground state, that may be easily prepared by optical 
pumping with 7r-polarizcd light. GHZ correlations are 
detected by a coincidence measurement of the photons 
emitted from the cavity and information on the atomic 
state obtained from a "Stern-Gerlach" device. Our re- 
sults are summarized in Sec. VIII. 

II. MODEL OF THE ATOM-CAVITY SYSTEM 

The experimental configuration we consider is similar 
to the one presented in [20] and is depicted in Fig. 1. A 
moving atom (transition frequency u a ) first encounters 
the waist of a quantized cavity field and subsequently 
a pump-laser field injected perpendicular to both cav- 
ity axis and the atom's direction of propagation. Cavity 
and pump fields overlap partially, so that there is a re- 
gion where in the rest frame of the atom the pump field 
increases as the cavity field decreases. The cavity field 
is described by two transverse-polarization states of a 
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FIG. 1. Setup of the two-polarization adiabatic passage 
scheme. With spatial quantization along the z axis, the cav- 
ity supports a a + - and a degenerate cr~-polarized mode. The 
atomic beam successively traverses the cavity waist and the 
7r-polarized pump laser beam, which is oriented perpendicular 
to the direction of the atoms and the cavity axis. 



single longitudinal mode. With the cavity axis as the 
direction for spatial quantization, the polarizations are 
taken to be er + and a~ and the resonance frequencies 
of the corresponding modes denoted as u) c ±. The abil- 
ity to accomodate arbitrary polarization of the cavity 
field distinguishes the model presented here from that of 
Refs. [20,23]. It not only allows us to investigate polar- 
ization entangled states of the cavity field, such as the 
maximally entangled GHZ states, but also leads to an 
improved scheme for the generation of single-polarization 
Fock states in a realistic experimental situation. 

We assume here that the atomic levels involved are 
the Zeeman substates of ground and excited levels F g 
and F e , although more general configurations can also 
be employed [22,25,26]. In most of the paper we con- 
centrate on the case F g — F e , where the transition 
(m g = 0) — ► (m e = 0) is electric dipole forbidden, so that 
a 7r-polarized pump beam does not couple to atoms in the 
Trig — state of the lower Zeeman manifold. Therefore, 
if by a suitable interaction with pump field and cavity 
mode, atomic population is pushed towards levels with 
decreasing \m g \, the process will terminate once the level 
(m g = 0) is reached, independent of the initial atomic 
distribution. Note that the coupling coefficients among 
the various Zeeman states for the F g <-> F e transition are 
symmetric under a change of sign of the magnetic quan- 
tum number, which is important for the generation of 
states involving photons of different polarization. 

The maximum number of photons that can be gener- 
ated with the scheme is achieved when the atom is pre- 
pared in the outer Zeeman ground states (m g — ±F g ) 
or a coherent superposition of both. As shown below, 
ideal adiabatic passage in this case will generate F g (en- 
tangled) photons distributed among both cavity modes. 
With cesium (nuclear spin I = 7/2, F g = 3 or 4) as an 
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FIG. 2. Scheme of the (F g = 3) and (F e = 3) hyperfine 
levels, showing the transitions relevant for adiabatic passage 
starting from two different initial preparations of Zeeman sub- 



levels (indicated by the solid circles): (a) 



-3; (b) su- 



perposition of m s 



-3 and m„ 



+3. The wide arrows 



symbolize the 7r-polarized (classical) pump transitions. The 
predominant cavity-induced transition from a given excited 
state e m [larger Clebsch-Gordan coefficient as indicated by 
the numbers in diagram (b)] is drawn as a solid, the weaker 
emission into the orthogonal polarization mode as a dotted 
arrow. Note that at the end of the passage the atom is al- 
ways found in the m 9 = level (open circle), even in the 
presence of two cavity polarizations. 



example, states with three or four entangled photons are 
therefore possible. Figure 2 shows the relevant pump 
and cavity-induced transitions in the case F g = F e = 3 
for two possible initial preparations of the Zeeman sub- 
levels. 

In our model, the cavity field is described by annihila- 
tion and creation operators a+,a' + for the c + -polarized 
mode and a_ , a[_ for er_ polarization, while the pump 
field is treated classically (i.e., as a coherent state in a 
regime of weak coupling). The atomic degrees of free- 
dom are described by operators A a , which arc composed 
of a sum over spin-1/2 operators for the Zeeman levels 
\F g ,m g ) and \F e ,m e ) with weights corresponding to the 
Clebsch-Gordan coefficients for dipole coupling with po- 
larization er, namely 

A a = ^2 \F g m g ){F g m g ;la\F e m e ){F e m e \, <r = 0, ±1. 



(1) 



As the atom moves across the cavity, it undergoes a time- 
dependent interaction with the fields, which is expressed 
by a variable coupling constant g(t) for the two quantized 
cavity modes and a Rabi frequency fi(i) for the classical 
pump field. The time scale is determined by the atomic 
velocity. Alternative time-dependent coupling schemes 
for the case of a stationary atom located in the cavity 
field [22,25,26] may be treated with the same formalism. 
With the above definitions, the Hamiltonian for the 
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atom and the two fields in a frame rotating at the fre- 
quency loi of the classical field is given by 

Hi nt (t) = H5 + a\_a + + hS-a^_a- 

-ihg(t) (alA +1 - A^a+^j 

+iMl(t) (Ao - 4) . (2) 

Here, S± = lu c ± — ui is the detuning between the cavity 
frequency iv c ± for the a± -polarized mode and the classi- 
cal field, which is assumed to be at resonance with the 
atomic transition {uo a = wj), throughout the present anal- 
ysis. 

A basis for the Hilbert space of the atom-cavity sys- 
tem is provided by tensor products of the atomic states 
and the number states for each polarization of the cavity 
mode. We use the notation 

\x m , n+, rc_) = \x m ) ® \n+) ® |n_), (3) 

with x = e, g specifying excited and ground state, m = 
—F x , ...,F X the number of the magnetic sublevel, and 
n± the photon number in the respective cavity mode. 

In this paper the initial preparation of the system is 
not described explicitly. We assume that before an atom 
enters the interaction zone both modes of the cavity field 
arc in the vacuum state and the atom is prepared in an ar- 
bitrary coherent superposition of Zccman ground states. 
This may be achieved by radio-frequency pumping [27] 
or a series of optical Raman pulses [28] . 

In our initial discussion in the following section, atomic 
spontaneous emission and cavity decay will be neglected. 
They will be included in Sees. IV and V in a master 
equation approach using the quantum Monte Carlo wave- 
function technique. 

III. DRESSED STATES AND DARK SPACE 

Fundamental aspects of the system evolution in the 
two-polarization adiabatic passage scheme may be in- 
ferred from the structure of the eigenvalues and eigen- 
vectors of Hi n t (t) during the transit of a single atom. As 
a specific example we consider the (&S\/2,F g = 3) — > 
(6P1/2, F e — 3) transition in cesium (cf. Fig. 2) at 852.36 
nm. For the system parameters we have chosen values 
reached in recent cavity QED experiments with high fi- 
nesse optical resonators [29-33]. In particular, a peak 
atom-cavity coupling g — 25T on the resonator axis was 
used with r/27r = 5 MHz being the spontaneous emis- 
sion rate for the cesium D2 line (i.e., T^ 1 is the popula- 
tion lifetime). Both cavity modes and the pump beam 
are assumed to have a Gaussian transverse shape with 
a full width at half maximum (FWHM) w = lOT' 1 , 
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FIG. 3. (a) Time evolution of the coupling g(t) and Rabi 
frequency Q(t) of the pump field, which are assumed to be 
Gaussian pulses with FWHM 10r~ ; they are centered at 
t — 17r _1 and t = 23T~ , respectively, and have amplitudes 
go — 25F and fio = 50T; (b) energy eigenvalues of _ffi nt in 
the single-polarization case; (c) energy eigenvalues for two 
polarizations. Atom, cavity modes, and pump field are at 
resonance (ui a — lj c+ = uj c - = u>i). 

but the pump beam has its center displaced along the 
atomic beam by an amount of 0.6u> relative to the cav- 
ity mode [Fig. 3(a)]. The maximum Rabi frequency fio of 
the pump field is chosen as twice the coupling constant go 
(i.e., f2o = 50r for the above specified coupling). These 
values were found to yield optimum results in the single 
polarization case [20,23] as well as in the present work. 
Small variations in the shape and intensity of the field 
distributions on the order of a few percent, however, do 
not lead to substantially different results. 

With the two cavity modes in the vacuum state and 
the atom prepared in the lowest Zeeman ground state 
(m g = —3), the initial wave function in the notation of 
Eq. (3) is 

|fi„i> = |<7-3,0,0). (4) 

The subsequent analysis will be restricted to the manifold 
of states that are coupled to |*i n i) by the interaction 



4 



Hamiltonian Hi nt . 

In order to compare the system behavior for the cases 
of one and two polarizations, we first calculate the eigen- 
values of Hi nt when only the tr_ mode interacts with 
the atom, corresponding to the single polarization anal- 
ysis of Refs. [20,23]. To this end the terms involving 
a+-operators are omitted from the Hamiltonian (2). Fig- 
ure 3(b) shows the eigenvalues obtained as a function of 
transit time for typical experimental parameters and a 
cavity detuning 5- = 0. There are seven nondegenerate 
energy levels that are also referred to as dressed states 
since they describe the energy of atom and the surround- 
ing fields. The dressed state employed for the adiabatic 
passage is the one that does not couple to excited atomic 
states and hence does not experience an ac-Stark shift 
due to either the cavity or pump fields (dark state). Its 
energy eigenvalue is E = and the corresponding eigen- 
vector in the basis (3), up to a normalization factor, is 

\Eo) = -Vr% 3 (t)| 3 _3, °. °> + V45.9 2 0, 1) 

-Vl8g(t)n 2 (t)\g- U 0,2)+n 3 (t)\g ,0,3). (5) 

This state has the property that if the system is pre- 
pared in |<7_3,0,0) in a region of f2 w 0, then slowly 
increasing Q(t) while reducing g(t) to zero (g/£l ~ 0) 
will adiabatically transform the initial state to \go,0,3), 
i.e., a state with three cr_ photons in the cavity mode. 
As is apparent from Eq. (5), the creation of each photon 
is accompanied by an increase in the magnetic quantum 
number m g of the atom. Since for a 7r-polarized classi- 
cal laser field no more pumping can occur once the state 
m g = is reached, exactly three photons will be found in 
the cavity mode at the end of the adiabatic passage. This 
method has been analyzed as a scheme for the generation 
of Fock states and for more general superposition states 
of the field arising from superpositions of initial Zccman 
states in the F g mainifold [20,23]. 

These conclusions arc no longer valid if a second, or- 
thogonally polarized cavity mode near resonance is taken 
into account, as a calculation of the eigenvalues of iJi nt for 
two cavity polarizations with equal frequencies demon- 
strates. According to Fig. 3(c), the resulting level struc- 
ture is considerably more complex than in the single po- 
larization case, with a large number of levels undergoing 
crossings and anticrossings. This is due to the fact that in 
contrast to the single-polarization case there is no intrin- 
sic limit to the number of photons that may be deposited 
in the cavity mode, as will be explained below. The nu- 
merical simulations (Sec. IV) show that not all the levels 
accessible to the system contribute equally to the dynam- 
ics. For example, a cutoff of six photons in both cavity 
modes was found to produce sufficiently accurate results 
and was used in Fig. 3(c). 

Another important implication of the second polariza- 
tion is that the dark state at the energy E = becomes 
highly degenerate, the degree of degeneracy being limited 
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FIG. 4. Fundamental processes of the dynamics within the 
dark space, (a) Scattering of a pump photon into the <r_ cav- 
ity mode; (b) scattering of a pump photon into the cr+-cavity 
mode; (c) polarization flip of a cavity photon. 

only by the cutoff photon number. There is no longer a 
single dark state but rather a dark space in which the 
dynamics evolve in the adiabatic case. Consequently, no 
simple final state of the atom-cavity system is expected 
to be generated by the adiabatic passage, but rather a 
complex superposition state in the dark subspace. 

In the single polarization case only one fundamental 
process is involved in the evolution of the dark state (5), 
in which a photon from the pump field is scattered into 
the <r_ cavity mode, increasing the magnetic quantum 
number by one [Fig. 4(a)]. The process repeats until the 
atom reaches the Zeeman ground state with m g = 0. 
As will be analyzed in more detail below, with two po- 
larizations present the final state of the adiabatic pas- 
sage still has the atom in m g = 0. However, two addi- 
tional processes contribute to the dynamics within the 
dark eigenspace of energy E = 0. In one process, a 
linearly polarized photon from the pump field is scat- 
tered into the cr + -polarized cavity mode, accompanied by 
a decrease in the magnetic quantum number of the atom 
[Fig. 4(b)]. As a secondary consequence, the system has 
to go through an additional cycle of the type shown in 
Fig. 4(a) to reach the final atomic ground state m g = 0. 
Effectively the system has gone through a loop in the 
atomic level scheme, creating a pair of a + /er_ photons in 
the process. Depending on the system parameters, sev- 
eral of these loops may occur. Clearly, loop processes 
must be avoided if a pure number state is to be gener- 
ated, as they lead to superposition states of different (and 
somewhat uncontrollable) photon numbers in the cavity. 

The second additional process is the destruction of a 
cr+ photon followed by the creation of a ct_ photon (and 
conversely) with no pump photons involved [Fig. 4(c)]. In 
this case the photon number in the cavity is not changed, 
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but the polarization of one photon is flipped. At the 
same time the atom changes its Zeeman substatc by two. 
Due to this process, states beyond the m g = Zeeman 
ground state may become temporarily occupied, which is 
impossible in the single polarization case. The branch- 
ing ratio determining the relative contribution of each 
of the three processes of Fig. 4 (at resonance) depends 
on three factors: the occupation of the cavity mode (in- 
fluencing induced emission), the Clebsch-Gordan coeffi- 
cients of the relevant transitions (cf. Fig. 2), and the tem- 
poral structure of the couplings and fields. By prepar- 
ing the atom in the lowest Zeeman state it is guaran- 
teed that the first transition is of type (a), enhancing 
the probability for further emissions with the desired cr_ 
polarization. The Clebsch-Gordan coefficients also favor 
transitions towards decreasing \m g \ so that the adiabatic 
passage method with two polarizations is still expected 
to work, albeit with an efficiency degraded by the pro- 
cesses of Figs. 4(b,c). To obtain pure Fock states of the 
cavity field for a given initial Zeeman substate, a method 
to suppress these undesired transitions must be found. 
We will demonstrate in the following that a finite detun- 
ing of the cavity from resonance effectively achieves this 
end. 

To analyze the structure of the dark space manifold, 
it is convenient to parametrize the vectors spanning the 
E = cigenspace by the number of photons contained. 
As photons arc transferred to the cavity during adiabatic 
passage, the initial occupation of the modes (when the 
pump field £1 = and the atom is in the state g_ 3 ) 
is used to characterize each eigenstate unambiguously. 
Only such states with an even initial photon number con- 
tribute, corresponding to the fact that nonadiabatic cou- 
pling to states with higher photon number only occurs by 
the loop process described above (which adds two pho- 
tons to the system). For example, for the initial atomic 
state | <7 3) , explicit expressions for the eigenstates corre- 
sponding to total initial photon numbers zero, two, and 
four are given by 



45< 7 2 (i)Q(*)|< 7 _ 2 ,0,l) 



Eo) = N [-V15ff d (i)|«7_3,0,0> 



\Ei) 



'I8 5 (t)fi 2 (t)| 5 -i, 0, 2) + n 3 (t)\g , 0, 3)] (6a) 
= M [\/60<7 3 (*)l<7-3,l,l) 
-vW(W)l</-2,l,2) 
+V24g(t)n 2 (t)\g^, 1, 3) - il 3 (t)| 3 o, 1,4) 
+\/l8< ? 3 (i)|< ? _i, 0, 2) - V36g 2 (t)n(t)\g 0l 0, 3) 
+V6g(t)Q 2 (t)\g +1 ,0A)] (6b) 
-M 2 [-\/l50 5 3 (t)| 3 - 3 ,2,2) 
+V^50ff 2 (t)fi(i)|<7- 2 ,2,3) 
-V36g(t)n 2 (t)\g- U 2, 4) + f> 3 (t)| 3o , 2, 5) 
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FIG. 5. Time evolution of the first two eigenvectors with 

degenerate eigenvalues (E = 0) in terms of the basis state 
1 1 2 

coefficients d m , n+ ,n_ • They are defined by the expansion 

\Ek{t)) = E m ,„± d ^,W,™- Wl£™> n +> rt -> t cf - E q- ( 6 )]- ( a ) 

Eigenvector with no er + contribution, corresponding to the 
dark state of the single-polarization case; (b) eigenvector with 
one 0+ photon in the cavity; note the additional contributions 
obtained by polarization flips from a + to ct_ (thin lines). 



12 ff (t)fi 2 (t)| ff+ i,l,5> 
15 5 (t) 3 | 5+1 ,0,4) 
15.g(t) 2 O(t)| ff+2 ,0,5>] 



(6c) 



-vW(*)l«7-i, 1, 3) + V90g 2 (t)n(t)\g 0l 1,4) 



with the symbols A4 denoting normalization constants. 
Note that \Eq) is identical to the single polarization dark 
state of Eq. (5). 

The evolution of the first two eigenstates (|£b) and 
\E\}) in terms of the basis state components \g m , n+, n_) 
is plotted in Fig. 5. As the adiabatic passage proceeds, 
each state acquires exactly three additional er_ photons 
due to the process of Fig. 4(a). Moreover, occupation is 
redistributed by the polarization flips shown in Fig. 4(c). 
For a change of eigenstate to occur, the system must un- 
dergo a loop process involving the generation of a er + / er_ 
photon pair. These latter events are the reason for the 
failure of standard adiabatic passage to produce a pure 
Fock state in the two-polarization case with degenerate 
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frequencies. 

Possible final states of the system may be deduced from 
Eqs. (6) by taking the limit g(t) — > 0. As expected from 
the vanishing coupling of the pump field in the central 
Zeeman state, the atom emerges from the adiabatic pas- 
sage always in the m = ground state as in the single 
polarization case. This is an important fact, because it 
ensures that the atom is not entangled with the field. The 
cavity modes are found in a superposition of the states 
|0, 3), 1 1,4), 1 2, 5), . . ., so that the general expression for 
the final state of the system is 



| So) ® ^c fe |fc,3 + fc). 



(7) 



k=0 



Here, the notation |n+,n_) is used to denote number 
states of the two-mode cavity field. 

As a consequence of Eq. (7), along with the desired a_ 
Fock state |0, 3), there is a finite probability for additional 
pairs of cr+ and <r_ photons in the cavity field. The size of 
those contributions, i.e., the modulus of the coefficients 
Cfe in Eq. (7), must be determined from the dynamics 
generated by H mt (see Sec. IV). However, an estimate 
for |ci| 2 can be obtained by restricting the system to the 
subspace spanned by the two lowest eigenvectors of the 
dark space, with \E ) being the state initially occupied 
and \E{) the first state with a a + polarization contribu- 
tion. As the eigenvectors in Eqs. (6) are time dependent, 
diabatic Landau-Zener-type transitions between the as- 
sociated eigenstates occur. The transition probability is 
obtained by integrating the Schrodingcr equation for a 
two-level system with time-dependent basis states: 



iftj^ - hJ\ [co(t) \Eo(t)) + ci(t) |£i(i)>] = 0. (8) 



Starting in state \Eq) at t — ► — oo, we derive the probabil- 
ity | cx | 2 for a diabatic transition to state \E\) at t — > oo 
to be given by the following expression: 



Ci =1 sin 



[Ei \ -§t\E ) e -i[ El (T)-E {T)]T/n dT 



(9) 

From Eqs. (6) one finds that the matrix elements ap- 
pearing under the integral do not vanish, as the eigen- 
vectors and their derivatives are not orthogonal [because 
of the assumed classical field for 0(f)] . The size of the 
relevant scalar products as a function of time and hence 
the diabatic transition probability peak in the middle of 
the atomic transit. For the parameters of Fig. 3, there is 
a 22% probability for the system to undergo a Landau- 
Zener transition during the adiabatic passage and end 
up in the state |so, 1,4). The model of two degenerate 
dark states thus demonstrates that adiabatic passage as 



described in [20] is not well suited to generate number 
states in a cavity if the two polarizations of the cavity 
field are closely spaced in frequency. 

For the discussion of a method to suppress diabatic 
transitions to states with the wrong polarization, it is im- 
portant to note the dependence of the integral in Eq. (9) 
on the phase factor. When \E$) and \E\) are degenerate, 
this phase factor is obviously zero. There is, however, a 
simple way to lift this degeneracy: a small detuning 5 of 
the cavity with respect to the pump field and the atomic 
transition frequency (assumed to be equal for both polar- 
ization modes, i.e., S + = 5_ = S) will split the dark space 
levels into nondegeneratc eigenstates, according to their 
total photon number. Therefore, these new eigenstates 
are exactly the states given by Eq. (6). The correspond- 
ing dark space eigenvalues form a series of energy levels 
separated by approximately 26, since neighboring levels 
differ by two photons. Through adiabatic passage each 
state acquires exactly three photons, so that all eigenval- 
ues undergo an increase of 3<5 during the atom's transit 
through the cavity. This is apparent in the central part 
of Fig. 6. 

When the degeneracy of the dark space is thereby 
lifted, the phase factor apppearing in Eq. (9) no longer 
vanishes, but leads to oscillations of the function under 
the integral. As a consequence, the transition probabil- 
ity | ci | 2 will drop when 6 is increased. This simple model 
is in excellent agreement with the results of a numerical 
integration of the full interaction Hamiltonian (2). 

With the transition rate strongly suppressed at large 
5, a system prepared in the |s-3,0,0) state will adiabat- 
ically follow the \E ) dressed state in the dark space and 
end up in the |so, 0, 3) state, i.e., a Fock state of the cav- 
ity field with pure cr„ polarization. In fact, the structure 
of the phase factor in Eq. (9) implies a new adiabaticity 
condition for the energy splitting of the dark state mani- 
fold and the time T for which the interaction is effective: 



\Eh. - Eq\T rs 2\S\T ^ 1. 



(10) 



This requirement must be met in addition to the usual 
(single-polarization) conditions for adiabatic following 
[20,34] 



n T»i, 2so V / ™± + 17 1 » 1, 



(11) 



which minimize the probability for diabatic transitions 
to ac-shifted (nondark) dressed states. Since for typical 
parameters 6 <C Uoi9o, the new condition (10) is more 
stringent than the standard adiabatic condition (11), so 
that considerably slower atoms (by a factor of go / 5) are 
needed in the two-polarization case. 

In the discussion, so far it was assumed implicitly that 
the energies of nondark dressed states (ac-Stark shifted 
on the order of go and Oo) and the dark space eigenvalues 
(on a scale of \5\) are well separated, as is the case in the 
middle of the atomic transit. However, this assumption is 
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FIG. 6. Energy eigenvalues of H ln t as a function of transit 
time for <5 = 0.5I\ The other parameters are the same as 
in Fig. 3. The bold line represents the single-polarization 
dark state \Eo). While at the center the splitting of the dark 
space manifold is well separated from all other levels, the level 
structure at the edges is rather complex. Note the degeneracy 
at Ft « 34. 



not valid at the leading edge of the cavity mode and the 
trailing edge of the pump field, where the ac-Stark shifts 
of the energy levels become comparable to the splitting 
of the dark states (g(t),Q(t) « \S\). In this regime, a 
complex structure of the eigenvalues ensues, as is appar- 
ent from the edges of the dressed-state diagram of Fig. 6. 
An enlarged view of the energy levels for the atom enter- 
ing the cavity mode (Fig. 7) shows that the lower dark 
state \Eo), which is initially occupied, undergoes a large 
number of crossings with nondark levels. Most of these 
crossings are actually avoided, as illustrated in the blow 
up of the first crossing in Fig. 7. Since the synthesis of 
pure Fock states requires that the system follows the low- 
est dark state, transitions to any intersecting states must 
be avoided, i.e., all crossings must be passed diabolically. 
This condition puts an upper bound on the products ap- 
pearing in Eqs. (10), (11), limiting them to a finite range 
of values. 

A quantitative evaluation of the eigenvalue spectra 
shows that the effects of level crossings are different when 
the atom is mainly exposed to the cavity field (in the ini- 
tial phase of the atomic transit with g/Sl ^> 1) and when 
the interaction is primarily with the pump field (in the 
final phase with g/Q, <C 1). In the former case, it should 
be noted that only two-level crossings occur. The largest 
level splitting is found at the first crossing, which is the 
one shown in the lower part of Fig. 7. Corresponding 
to the small coupling at the edge of the cavity mode, 
it has an energy splitting orders of magnitude smaller 
than T (3.5 x 10~ 4 r in the example of Fig. 7). In order 
to pass this crossing adiabatically, the atom would have 
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FIG. 7. Close-up view of energy eigenvalues of Hint at 
the leading edge of the cavity mode for the parameters of 
Fig. 6. The enlargement of the first crossing of the lowest 
dark state (with state |e_2,0, 1)) shows that the crossing is 
actually avoided on a very fine scale. 



to have a velocity 10 6 times lower than that taken for 
Fig. 3. Therefore, for realistic atomic velocities the dia- 
batic condition is easily fulfilled on entry and the atom 
will pass the crossings as if the levels were not coupled. 
This result is confirmed by the numerical calculations of 
the following section. 

A different situation occurs at the trailing edge of the 
pump field. At a time t such that the classical Rabi 
frequency has decreased to fl(t) = 4|<5|, the ac-shifted 
dressed levels become degenerate with the series of dark 
states (neglecting coupling to the cavity mode) , as is ap- 
parent from Fig. 6 and shown in more detail in Fig. 8 for 
the lowest dark state. The phenomenon is found in the 
single- [Fig. 8(a)] and the two-polarization case [Fig. 8(b)] 
at finite detuning. The degeneracy is lifted due to the in- 
teraction with the cavity mode(s), so that a multiple-level 
avoided crossing is created. For positive detuning 8 and a 
single polarization, the multiplicity is given by the num- 
ber of dressed states with positive ac-Stark shift [three 
levels crossing the dark state in the example of Fig. 8(a)]. 
In order to emerge from the crossing in the dark state, 
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FIG. 8. Dressed level structure at the trailing edge of 
the pump pulse where fl(t) ~ 4|<5|. A degeneracy of the 
atom-pump interaction leads to a multilevel avoided cross- 
ing, (a) Single-polarization case; (b) two-polarization case. 
The parameters are the same as in Fig. 6. 



the system would have to undergo diabatic jumps across 
three levels with a substantial energy splitting (on the 
order of 0.02r for 5 = 0.5T in the example of Fig. 8). 
Consequently, there is a finite probability for transitions 
to states outside the dark state manifold. These partly 
adiabatic transitions at the trailing edge of the pump 
field limit the detuning 5 that may be applied to achieve 
adiabatic conditions at the center of the interaction re- 
gion. They also set an upper limit on the transit time 
T and hence on the products of Eqs. (10) and (11), as 
pointed out earlier. 

For two polarizations, an additional degeneracy occurs 
due to contributions from states with higher photon num- 
ber. In this case dressed states with negative and positive 
Stark shift are involved in the crossing [Fig. 8(b)]. Note 
that a similar multiple-level crossing does not occur when 
the atom enters the cavity, as the dressed energies of the 
quantized cavity mode are nondegenerate due to their de- 
pendence on the photon number. Hence the restrictions 
for detuning and time of flight are more stringent during 
interaction with the classical field in the final phase of 



the atomic transit. 

In the following sections, the preceding conclusions 
drawn from the eigenvalue spectrum of the Hamiltonian 
in Eq. (2) and a simple two-level model will be put on a 
quantitative basis by integrating the master equation of 
the system, including dissipative processes in the treat- 
ment. A suitable method is the quantum Monte Carlo 
wave-function technique. The results will be employed 
to describe a strategy and calculate the efficiency for the 
synthesis of pure |0,3) or |3, 0) states of the cavity field 
or superpositions of both. 

IV. QUANTUM MONTE CARLO SIMULATIONS 

In the previous analysis only the coherent interaction 
of the atom and the two cavity modes was considered. A 
realistic model of the system also has to take into account 
dissipation that enters due to spontaneous decay of the 
atom and damping of the cavity mode. The dissipative 
atom-cavity system is described by the following master 
equation for the density operator p(t) of the system: 

dp i 

— = -- [H cS , p] + 2K(a + pa\_ + a_pat) 

+r J2 a °p a I ( 12 ) 

cr = 0,±l 

Here, i? e ff is a non-Hermitian effective Hamiltonian in- 
cluding decay terms given by 

H eS = H int - iK(a\a + + aLa-) - i- ^ A\A a . (13) 

(T=0,±1 

Solutions of the master equation may be obtained by 
integrating Eq. (12) numerically. However, the size of 
the Hilbert space accessible to the system in the two- 
polarization case is rather large, even if its basis is trun- 
cated at moderate photon numbers. For example, for 
atomic levels with F g = F e = 3 and a cutoff photon 
number of seven in both cavity modes, the dimension N 
of the Hilbert space is2x(2x3+l)x8x8 = 896. 
The computational expense involved in density matrix 
calculations, which scales as A^ 4 , makes integrating the 
master equation impractical. A more suitable approach 
is to simulate the master equation by using Monte Carlo 
wave- functions [35-37] . Here, only the Schrodingcr equa- 
tion with the non-Hermitian Hamiltonian H c s has to be 
integrated and the problem size scales as A^ 2 . Dissipa- 
tion is taken into account by quantum jumps, which oc- 
cur at random times with a probability determined by 
the decay rate of the wave-function norm. The jumps 
are implemented by applying so-called collapse operators 
to the system wave- function and normalizing the result. 
Between subsequent jumps, the system evolves accord- 
ing to the Schrodingcr equation with the non-Hermitian 



9 



10 20 30 40 

Time t (in units of r 1 ) 

FIG. 9. System evolution of an atom and two degen- 
erate polarization modes, obtained from a Monte Carlo 
wave-function simulation including atomic spontaneous emis- 
sion. Shown is the occupation probability of the basis states 
of the system. The amplitudes of the Gaussian pulses are 
go = 25F (at t = 17F- 1 ) and Q = 50r (at t = 23r- x ), with 
a FWHM of 10r _1 . The cavity is at resonance with the pump 
field (S± — lo c ± — uii — 0). 

effective Hamiltonian H c ff. In this way, one obtains a 
wave function conditioned on a particular decay history. 
The system density matrix and operator expectation val- 
ues are obtained by repeating the calculation many times 
and averaging over the resulting trajectories. Details of 
the procedure are found in Refs. [35-38]. 

In addition to the potential for more efficient use of 
computational resources, the method offers the benefit 
that for each of the simulated trajectories a record of the 
quantum jumps that have occurred during the system 
evolution is obtained. In the case of Eq. (12), there are 
five possible collapse operators and hence five types of 
jumps given by 

V2na±, VTA,,, a = 0,±1, (14) 

corresponding to decay of the two orthogonally polarized 
cavity modes and spontaneous emission of a photon with 
polarization a, respectively. Decay of the cavity modes 
may be associated with photodetection events in an indi- 
vidual experimental run. Note that, when the polariza- 
tion of the photons leaking from the cavity is analyzed in 
a basis different from that of the intracavity field (as set 
by the spatial quantization of the atom), linear combina- 
tions of the field decay operators a± must be employed 
[39], which will be used in Sec. VI to analyze GHZ cor- 
relations between cavity photons. 

In a first set of simulations, only atomic decay is in- 
cluded, while cavity decay is still neglected. Ideally, as 
long as the evolution of the system is restricted to the 
dark space, the influence of spontaneous emission on the 
dynamics should be weak, since the upper atomic levels 
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Time t (in units of r ) 

FIG. 10. System evolution of atom and two polarization 
modes, at a finite cavity detuning 8 = 0.6r, showing nearly 
ideal adiabatic passage. The amplitudes and timing of the 
applied pulses are the same as in Fig. 9. 

are not excited. Hence, the conclusions of Sec. Ill should 
apply to the numerical results in this regime. 

We start by studying the system evolution for an 
F g = 3 — > F e = 3 transition in the case of atom, cavity, 
and laser fields on resonance (S = 0), i.e., the case of a de- 
generate dark space. The initial state is again |g_ 3 , 0, 0). 
Figure 9 shows the system evolution during the passage of 
an atom through the time-dependent cavity coupling g(t) 
and the laser field in terms of the basis state occu- 
pation. The probability maximum is shifted in turn from 
the initial state to [5-2,0, 1) to |g_i,0,2), before finally 
reaching \go, 0, 3). These are exactly the levels comprised 
by the lowest dressed state \E ) in the dark space man- 
ifold of Eq. (6), which should be the only ones occupied 
if the interaction was entirely adiabatic. Approximately 
halfway through the transit, however, contributions from 
levels |<7o,l,4) and \go,2,5) start to appear. These are 
due to the diabatic transitions to states with higher pho- 
ton number predicted from the two-level model analyzed 
in Sec. III. At the end of the transit, the atom is found in 
state \go) and the cavity is in a superposition of the states 
|0, 3), 1 1, 4), and [2, 5), etc., instead of the desired pure 
Fock state |0, 3), confirming the earlier result of Eq. (7). 

In the two-level model of Eqs. (8), (9) in Sec. Ill, a 
finite cavity detuning <5 was shown to lift the degeneracy 
of the dark space and lead to a strong suppression of di- 
abatic transitions. What remains to be investigated is 
the question of whether this mechanism is still effective, 
if the full system dynamics including atomic decay and 
level crossings in the entrance and exit phase of the inter- 
action are taken into account. In Fig. 10, a Monte Carlo 
calculation of the system evolution for a cavity detuning 
S± = 0.6T is shown. The pump field is still assumed to 
be on resonance with the atomic frequency (wj = u a ). 
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FIG. 11. Final state probabilities as a function of cavity 
detuning 8 + — <5_ — 8, showing the transition to an adiabatic 
regime for 8 > 0.5r. The parameters correspond to those of 
Fig. 9. 

As can be seen from the resulting final state, the finite 
cavity detuning leads to the desired strong suppression 
of states with ct+ photons and an amplitude for the ct_ 
Fock state, which is close to unity (99%). The atom re- 
mains in the lowest branch of the dark state manifold and 
is adiabatically transferred from |g_ 3 ,0,0) to \go,0,3)- 
Therefore, passage through a detuned cavity is adiabatic 
and well suited to generate a Fock state of one polariza- 
tion of the cavity field, even for a cavity supporting two 
degenerate polarization modes. By choosing other initial 
states (|<7_ m , 0,0), m = 1, 2), states with photon number 
m < F g may be obtained as well. Coherent superpo- 
sitions of Zeeman ground states with < m < F g are 
also possible, leading to superpositions of photon num- 
ber states. Thus an operationally simple modification of 
the work of Ref. [20] suffices to make a workable scheme 
for the generation of "arbitrary" field states for a single 
polarization mode. 

To study systematically the effect of cavity detun- 
ing 5 (assumed to be equal for both polarizations, i.e., 
5+ = 5- = S) on the efficiency of our scheme, we have 
calculated the final state as a function of 5, with the re- 
sult displayed in Fig. 11. As S is increased from zero, 
the probability for contributions from basis states with 
the "wrong" cavity polarization (cr + ) drops to zero, while 
the pure cr_ state becomes occupied with 99% efficiency 
around S ~ 0.6r. 

In contrast to the results of the two-level model, the 
probability for a pure cr_ Fock state drops again slightly 
when 5 is increased to values larger than 0.6r. The rea- 
son is the dynamics in the initial and the final phase of 
the atomic transit (for example, at the multilevel cross- 
ing shown in Fig. 8), which becomes partly adiabatic 
if S is chosen too large. As a consequence, population 
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Cavity detuning 5 (in units of r) 

FIG. 12. Probability for a final state |<7o,0, 3} as a func- 
tion of cavity detuning 8 for three coupling constants in the 
case of two-polarization adiabatic passage (solid lines). The 
single-polarization (a_) results are included as well (dotted 
lines). The timing of the pulses is the same as in Fig. 9 and 
Qo = 2go. Note that, for this figure, the full master equation 
was simulated in order to obtain an accuracy well below 1%. 



is transferred to levels outside the dark state manifold, 
leading not only to a finite probability for a decreased 
photon number in the cavity, but also to entanglement 
between the atom and the cavity field. This is illustrated 
in Fig. 11 by the rise of contributions from levels \gi~, 0, 2) 
(k = —3, . . . , 0) as 6 is increased beyond 0.6r. Note 
that this limitation of adiabatic passage also occurs in 
the single-polarization case and is not connected to the 
structure of the dark space. 

A blown-up plot of the probability for reaching the fi- 
nal state |<7n,0, 3) vs detuning 5 is presented in Fig. 12. 
The results for three different values of the coupling 
constant go show that the efficiency of the method ap- 
proaches unity when the coupling is increased (99.6% at 
go = 50T). The optimum detuning is 5 w 0.6T and de- 
pends on the coupling only weakly. For comparison, the 
single-polarization probability for a three-photon Fock 
state is also included in Fig. 12. The large deviation be- 
tween the two sets of curves for values below S = 0.QT 
clearly shows the inadequacy of the single-polarization 
model. For large detunings, on the other hand, the re- 
sults coincide, demonstrating that in this regime the two- 
polarization dynamics is well described by the adiabatic 
following of a single dark state. The common asymptotic 
behavior of the curves is another hint to the polarization 
independent mechanism populating nondark states at the 
multilevel avoided crossing occurring for large S (Fig. 8). 
Figure 12 suggests that, for an optimum value of S, the 
purity of the final state is only limited by the atom-cavity 
coupling go that may be achieved experimentally, if the 
interaction time is held constant. 
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FIG. 13. Probability for a final state \go, 0, 3} as a function 
of the detuning 8+ of the <r+ polarized cavity mode. The <t_ 
mode is assumed to be on resonance and a coupling go = 25r 
was used for both modes. The timing of the pulses is the same 
as in Fig. 9. 



In the calculations presented so far, the detunings of 
the two polarization modes of the cavity were always as- 
sumed to be equal (<5+ = 5-). In order to probe the 
sensitivity of the system to a finite mode splitting, which 
may occur due to cavity imperfections, we have varied 
the detuning 5+ of the er + -polarizcd mode while keep- 
ing the (7_ -polarized mode at resonance with the atom 
(5- = 0). Again, for the initial state |g_3,0,0), Fig. 13 
shows the probability to reach a hnal state \go, 0, 3) as a 
function of S + . It exhibits a gradual transition from the 
diabatic losses in the degenerate two-polarization case to 
the perfect adiabatic passage obtained for a single (cr_) 
polarization, as S+ is increased from zero. The single 
mode description applies when the a+ mode is far enough 
detuned from resonance for its influence on the system 
dynamics to be neglected, which is the case for detun- 
ings <5 + > T. Figure 13 suggests that, for the model 
of degenerate modes to be valid, the splitting should be 
small compared with the natural linewidth T. As long 
as S + < O.ir, the final state probability deviates from 
the degenerate case by less than a percent. Note that for 
realistic cavities a well-defined splitting in the a± basis 
may be difficult to achieve, because of a lack of control 
over the polarization basis in which birefringence of the 
mirrors occurs. 

At fixed coupling g , the efficiency of the state syn- 
thesis scheme may be increased by extending the time 
T during which the atom interacts with the cavity mode 
and the pump field, i.e., using slower atoms for the adi- 
abatic passage. In Sec. Ill it was found that adiabatic 
transitions to nondark states place an upper limit on T. 
There is, however, another restriction for the interaction 
time, which is set by the decay time of the cavity modes. 



This will be demonstrated in the following section by in- 
cluding cavity decay in the master equation. 



V. PREPARATION OF A 
GREENBERGER-HORNE-ZEILINGER STATE 

Adiabatic passage with a detuned cavity and an atom 
prepared in the \gs) level was shown in Sec. IV to syn- 
thesize Fock states of er_ polarization and effectively sup- 
press excitation of the <7+-polarized cavity mode. An 
important property of the proposed scheme is that it is 
entirely symmetric under a change of sign in the mag- 
netic quantum number and a simultaneous exchange of 
<r_ and a + polarization, as is apparent from the mas- 
ter equation (12) and the relevant Clebsch-Gordan coef- 
ficients. If, for example, the system is prepared in the 
state |g+3,0, 0), a Fock state with three <x + photons in 
the cavity (|3, 0)) will be generated. Another notable 
feature of the two polarization adiabatic passage for an 
F —* F transition is that at the end of the interaction the 
atom is always in the m g — Zeeman ground state (|<7o))i 
independent of the initially occupied Zeeman substates. 
Therefore, the atomic part of the wave function is always 
disentangled from the state of the cavity field and can be 
disregarded without introducing projection noise to the 
measurement of the cavity field. The method is not re- 
stricted to single Zeeman levels as initial atomic states. 
The atom may be prepared in a coherent superposition 
of Zeeman ground states, leading to superposition states 
of the cavity field as described for the single-polarization 
model in Refs. [20,23]. A particularly intriguing possi- 
bility in the two-polarization case is to use coherent su- 
perpositions of Zeeman sublcvcls with different sign of 
the magnetic quantum number m. In this case, quantum 
superpositions of photon states with different polariza- 
tion are synthesized, i.e., the adiabatic passage maps the 
Zeeman coherence of the atomic state to polarization en- 
tanglement of the two cavity modes. 

An important example of such field states are the num- 
ber states of n photons maximally entangled in polariza- 
tion, which are obtained for a superposition state of n 
photons in the <t + mode and n photons in the er_ mode. 
With the notation previously introduced for a two-mode 
cavity field, these states may be written as 
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(|n,0) + |0,n». 



(15) 



For the state |4 , 2), there is a close correspondence with 
the correlated two-particle states used to test Bell's in- 
equalities [1,40] obtained from local hidden- variable the- 
ories. Numerous experiments have shown these inequal- 
ities to be violated [41]. Yet even more interesting prop- 
erties emerge from the entangled state of three par- 
ticles, corresponding to a so-called GHZ state, as was 
introduced by Greenberger, Home, and Zeilinger [3]. Its 
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FIG. 14. Monte Carlo simulation of the preparation of a 
GHZ-type state in the absence of cavity decay. The am- 
plitudes of the cavity and pump field are go — 25T (at 
t = nF' 1 ) and fi = 50r (at t = 23F- 1 ) with a FWHM of 
10r _1 , the cavity detuning is S — 0.6F, and the initial prepa- 
ration of the atom is a coherent superposition of m g = 3 and 
m g — —3. (a) Occupation of the cavity modes (the distribu- 
tions are identical for the o+ and the cr_ mode); (b) atomic 
ground-state populations (the curves for +m g and — m g co- 
incide). The final state is a 99% pure state |GHZ). 



significance is based on the fact that, to test theories of 
local realism against quantum mechanics, a single set of 
observations of this state is sufficient "all or nothing"). 
This is in contrast to the statistical nature of the viola- 
tion of Bell's inequalities (for n = 2). 

There is a notable difference between the states (|n, 0) 
+ |0, n)) considered here and the original Bell-type states 
of the form (|+i, + 2 , • • • , +„) + hi, -2, • • • , -«))■ While 
in the former case all photons of a given polarization oc- 
cupy a single mode, in the latter case they are distributed 
over n distinct (i.e., distinguishable) modes (1,2,..., n). 
However, as will be discussed below, we restrict our at- 
tention to the subensemble of photoelectric detection 
events that produce one and only one detection event 
("click") at each of n sets of detectors. Hence, from 
the perspective of a post diction, the state |*„) for the 
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FIG. 15. Monte Carlo simulation of adiabatic passage for 
the same parameters as in Fig. 14, but with a cavity decay 
rate k — 0.2F. (a) Occupation of the cavity modes; (b) atomic 
ground-state populations. 



subensemble of distinct "temporal" modes generates cor- 
relations identical to those of the Bell-type states. More 
generally, by tailoring the interaction profiles g(t) and 
Q(t) of cavity and pump field, an improved timing of 
cavity decay events may be achieved, so that in fact suc- 
cessive quanta 1,2, ... ,n would emerge from the cavity in 
well defined intervals and would be routed in a predefined 
fashion to respective detectors, making the requirement 
of separate detection events even less restrictive. 

In the following, we give a detailed description of how 
the two-polarization adiabatic passage scheme may be 
extended to generate a GHZ-type state l^) (denoted 
as | GHZ)) with high efficiency. Towards this end the 
F g = 3 — > F e = 3 transition used for the generation of 
the single-polarization Fock state will be employed again. 
Note that states |\I>„) with n > 3 may be synthesized as 
well, if an atomic transition with a high enough angular 
momentum were to be used. For the F g = 4 — ► F e = 
4 hyperfine transition in cesium, for example, the state 
1*4) may be generated in exactly the same way as 1^3). 

The starting point for the generation of a GHZ-type 
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state is to create a coherent superposition of the two 
outer Zeeman ground states (m g — ±3). This task may 
be achieved by using suitably timed microwave [27] or 
Raman pulses [28]. Ideal adiabatic passage of an atom 
through the cavity modes with coupling g(t) and the 
pump field either by physical motion or by suit- 

able modulation of the coupling for a stationary atom in 
the cavity, would then generate the following mapping of 
ground-state Zeeman coherences to the cavity field: 

(Iff-s) + 1.9+3)) ® |0, 0) =► \go) ® (|0, 3) + |3, 0)) 

~ | 30 )®|GHZ) (16) 

and would thus leave the cavity in a maximally entangled 
state of the GHZ-type. The transitions relevant for this 
synthesis scheme are the ones shown in Fig. 2(b). 

We have tested the efficiency of such a GHZ-state 
preparation procedure numerically by simulating the rel- 
evant master equation using Monte Carlo wave- functions, 
again neglecting cavity decay in the first instance but 
with atomic decay included. The same parameters as in 
Sees. Ill and IV are used with a detuning 5 = 0.6T to 
achieve optimum purity of the final state. The resulting 
averaged time evolution is shown in Fig. 14. In Fig. 14(a) 
the occupation of the cavity mode states is plotted, while 
Fig. 14(b) shows the corresponding atomic ground-state 
populations. As expected due to the symmetry of the 
system, the time-dependent photon distribution of the 
two cavity modes is identical. The same is true for the 
occupation of Zeeman ground-state levels with the same 
modulus \m g \. At the end of the adiabatic passage, there 
is almost perfect polarization entanglement with a 99% 
probability for the cavity to be found in the GHZ state 
( 1 0, 3} + 1 3, 0) ) and the atom in state \go). In correspon- 
dence with the discussion of Sect. Ill, contributions of 
higher terms in the dark space manifold, which are again 
of the form (\k, 3 + k) + |3 + k, k}), k = 1,2,..., are 
strongly suppressed. As mentioned above, the efficiency 
could be further increased by using an even stronger cou- 
pling go or slower atoms. 

Beyond these encouraging results, a realistic simula- 
tion of our experiment has to take into account dissipa- 
tion of the cavity field, which has been neglected so far. 
In recent experiments cavity decay rates k/2tt ranging 
from 0.6 to 35 MHz have been achieved [29,31-33]. As a 
point of reference, here we take initially k/2it = 1 MHz, 
resulting in k/T = 0.2, although we will subsequently 
investigate the dependence of the state generation on k. 
The cavity decay time 1/k is then comparable to the adi- 
abatic passage time scale of maximum efficiency for state 
generation [cf. Fig. 3(a) for the time dependence of g(t) 
and fl(t)}, so that cavity decay will play an important 
role in the present scheme and must be considered in the 
simulation. 

We have repeated the Monto Carlo wave-function sim- 
ulation of the GHZ state generation of Fig. 14 using 



the full master equation (12), including the cavity de- 
cay terms with k — 0.2T. In Fig. 15, the resulting aver- 
age time evolution of the cavity mode occupations and 
the atomic ground-state populations is displayed. As ex- 
pected, cavity decay strongly affects the dynamics of the 
field modes. During the adiabatic passage, the probabil- 
ity for the excitation of photon states with n > is sub- 
stantially reduced compared to the case without cavity 
decay. The probability for one photon in a cavity mode 
never exceeds 25%. When the atom leaves the cavity, 
both modes of the cavity field have already decayed to 
the n = ground state at the rate k _1 [Fig. 15(a)]. Ob- 
viously, there is no time when the cavity modes approach 
a pure GHZ-type state. This result is in correspondence 
with previous studies of the generation of quantum su- 
perpositions in a lossy single-polarization cavity [20,23]. 

Perhaps surprisingly, Fig. 15(b) shows a completely 
different behavior for the atomic population. The occu- 
pation of the Zeeman states is almost unaffected by the 
presence of cavity decay. The evolution proceeds essen- 
tially undisturbed in close correspondence to the perfect 
cavity case shown in Fig. 14(b). The fact that cavity 
dissipation does not affect the atomic dynamics during 
the adiabatic passage was also observed in the single- 
polarization calculations of Ref. [23]. It implies that 
atomic coherence is transferred without "damage" as the 
cavity field undergoes a quantum jump to a lower dark 
state manifold [42]. The preservation of atomic coher- 
ence during the adiabatic passage even for a lossy cavity 
is an essential requirement for the detection of GHZ cor- 
relations between the photons leaking from the cavity, as 
will be explained in the following section. 

VI. DETECTION OF GHZ CORRELATIONS 

In principle, there are two different approaches to de- 
tecting nonclassical states of the cavity field. First, the 
field inside the resonator may be probed directly. This 
may be achieved, for example, by sending a second atom 
through the cavity and reversing the adiabatic passage 
scheme to map the photon number distribution of the 
cavity field to coherences in the Zeeman substructure 
[20,23,43]. Measurements of the atomic substates would 
then provide information on the state of both field modes 
(in principle the "complete" state [43]). The required 
timing of the pump and probe atoms would certainly be 
difficult to achieve for an atomic beam (with Poissonian 
fluctuations). However, in the optical domain, the most 
severe obstacle to this approach is cavity decay, since, as 
shown above, dissipation of the cavity field prevents the 
modes from ever reaching a pure GHZ state, so that any 
attempt to observe a signature of that state would fail. 

This state of affairs suggests observing the photons 
leaking from the cavity mirrors during decay of the field 
as the method of choice for detecting GHZ correlations, 
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especially since it affords the opportunity for tests of lo- 
cal realism with the "flying" photons that emerge from 
the cavity. Correlations between these emitted photons 
could readily be measured by coincidence detection, with 
these measurements performed in spatially separated re- 
gions as required by the locality assumption. 

Before presenting a more detailed proposal for an ex- 
perimental setup to probe the decay photons, we will 
first discuss the expected correlations between photons 
emerging from a lossy cavity. For simplicity we start by 
assuming that the two cavity modes already are in the 
state | GHZ) given by Eq. (15) for n = 3. Afterwards we 
will examine how the correlations are affected if dissipa- 
tion is present during the state synthesis. 

Greenberger, Home, and Zcilingcr have developed 
their argument for three arbitrary particles in a maxi- 
mally entangled state. Therefore, their analysis can also 
be applied to three photons in the state |0, 3) + |3, 0). To 
obtain an operator with two eigenvalues ±1, as consid- 
ered by GHZ, we define an observable equivalent to the 
spin of an atomic system. It is easy to verify that the fol- 
lowing operators satisfy the spin commutation relations 

L + = a} + a-, L- = a}_a + , L z = 7j( a + a + — aLa_). 

(17) 

The operators L± perform a change of polarization on 
a cavity photon, which is analogous to a spin flip. The 
observable corresponding to the particle analyzer in the 
setup proposed by GHZ is then given by the superposi- 
tion of L + and L_ with a variable relative phase shift fa 
namely 

L(fa) = e^' 2 L+ + e- 4 ^ 2 L_ , (18) 

with the required eigenvalues ±1. To arrive at a measure- 
ment procedure for L{<j>) in the case of photons emitted 
from the cavity, we first rewrite 

L(fa = aUfaa x (fa - al(faa y (fa. (19) 

The operators a x ((j>) and a v {fa) are annihilation operators 
of linearly polarized photons in the x and y directions of 
a coordinate system rotated by an angle fa/4. They are 
related to the cavity mode operators by 

a x (<f>) = a + + e^ 4 a-), 

a y (t) = ^=(e-^a + -e^a_). (20) 

Therefore, according to Eq. (19), L((j>) can be determined 
by measuring the difference in the occupation number of 
two orthogonal linear polarization modes in a rotated 
frame. 

To be consistent with the GHZ gedanken experiment, 
we must make the additional assumption that each of 




FIG. 16. Idealized experimental setup for the detection of 
GHZ correlations between photons leaking from the output 
port of the cavity. Each of the three decay photons is sent 
by router R to a particular analyzer, consisting of a half-wave 
plate, a polarizing beamsplitter PBS, and a pair of detectors 
it D i5 counting photons of linear polarization £ (£=x,y) in the 
basis determined by the wave plate. 

the three photons from the state |GHZ) is delivered to a 
separate analyzer. In this case there are only two possible 
outcomes of a measurement of L(fa): a value of +1 is 
obtained if an x-polarized photon is detected, a value of 
— 1 for a y-polarized photon. 

Figure 16 shows an idealized experimental setup for the 
detection of GHZ correlations. We assume well-defined 
generation times for each of the three photons, as might 
be achieved in principle by a succession of suitable shifts 
applied to the atomic levels in concert with a somewhat 
revised time history for g(t) and Q(t) during the adia- 
batic passage. Each photon emitted from the cavity is 
then directed to one of three analyzers. In front of each 
analyzer (k — 1,2,3) there is a half-wave plate rotating 
the respective linear polarization basis by angle faj4. 
The polarization [x(<j>k) or y(4>k)\ is measured with a lin- 
ear polarizing beamsplitter (PBS) and photon counters 
at the two output ports. The analyzers are assumed to 
be positioned far from the cavity and from each other, 
so that the individual measurements may be regarded as 
causally separated. 

To arrive at the correlations discussed by Greenberger, 
Home, and Zeilinger, the product of the three detector 
results is considered. It is straightforward to calculate its 
expectation value if the photons are in the GHZ state 

(GHZ|Li(^i)L 2 (02)L 3 (^3)|GHZ) = cos (fa + fa + fa) , 

(21) 

where again the subscript identifies the analyzer. 

The most interesting situation occurs when the detec- 
tor angles are set such that (fa + fa + fa) is an integer 
multiple of tt. In this case the average value is either +1 
or —1, which coincides with the maximum and minimum 
possible value, respectively, since any individual product 
of detector outcomes is either +1 or —1. This implies 
that each individual measurement of triple correlations 
yields the same result. In this case the photon triplets 
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show perfect correlations in a sense that knowing the re- 
sults of two detectors allows one to predict the outcome 
of the third measurement with certainty. 

In Ref. [3] it was shown that the triple correlations of 
particles prepared in a GHZ state violate theories based 
on the assumption of local realism in a rather striking 
way. This may be demonstrated by considering the cor- 
relations for the following choice of detector angles: 

(GHZ| Li(f )L 2 (f )L 3 (0)|GHZ) = -1, (22a) 

(GHZ| Li(f )L 2 (0)L 3 (f )|GHZ) = -1, (22b) 

(GHZ| Li(0)L 2 (f )L 3 (f )|GHZ) = -1, (22c) 

(GHZ| Li(0)L 2 (0)£ 3 (0) |GHZ) = +1. (22d) 

If the result of each detector was independent of the an- 
gles 4>k of the other detectors (as required by the locality 
assumption for widely separated detectors), the expecta- 
tion values of Eqs. (22) would factor into the product of 
three functions Ak(<fik) that can take values ±1: 



(GHZ|L 1 (0i)i 2 (0 2 )L 3 (<^ 3 )|GHZ) 



(23) 



Multiplying the four left-hand sides of the set of Eqs. (22), 
a product expression is obtained in which every factor 
Ak((j)) (k — 1...3;4> — 0, 7r/2) appears exactly twice. 
Therefore, using A\ = 1), a result of 1 is expected for 
the overall left hand side product, if local realism holds. 
On the other hand, the right-hand sides of Eqs. (22), 
which were obtained quantum mechanically, multiply to 
give —1, clearly contradicting local realism. 

We now return to the situation encountered in our 
experiment for realistic parameters of coupling, transit 
time, and cavity decay time. It has been shown in Sect. V 
that, even for the best Fabry-Perot cavities available [44] 
(excluding for the moment the whispering gallery modes 
of quartz microspheres [45-47] ) , decay may occur well be- 
fore the adiabatic passage is terminated so that the intra- 
cavity field never reaches the state |GHZ). Further, since 
there is not a well-defined time history for the times of 
the photon emissions, the operator expectation value (21) 
is not necessarily a relevant quantity. However, correla- 
tions between the photons emitted during adiabatic pas- 
sage may still be probed. The important question then 
becomes whether or not they show a similar behavior as 
that described by Eq. (21). 

The Monte Carlo wave- function method offers a conve- 
nient way to evaluate correlations between photons emit- 
ted from the cavity. For each simulated trajectory a clas- 
sical record of the quantum jumps, which have occured 
during the adiabatic passage, is kept that contains in- 
formation on the polarization of the decay photon cor- 
responding to the jump. Averaging over many runs, 
photon-counting distributions and correlation functions 
for the detectors involved are readily calculated. To dis- 
tinguish between the six possible detection channels (two 



for each of the three analyzers), the collapse operators for 
the two cavity modes, which appear in Eq. (12), must be 
replaced by a new set of annihilation and creation oper- 
ators. As each decay of the cavity field can be associated 
with a click in a certain detector, six new collapse oper- 
ators are needed. They involve the operators defined in 
Eq. (20) in the rotated linear polarization bases: 



(24) 



2^a x (<t> k ) = ^{a+e-^ 4 + a_e^/ 4 ); 
2^a v ((f> k ) = i^(a+e- l<t "'/A - a_e^/ 4 ); 
k = 1,2,3. 

The new master equation to simulate is 



<7 = 0,±1 



+ y 22 { a x{4>k)pal{4>k) +a y (<j)k)pal(<t>k)) ■ (25) 

fe=l,2,3 

In the original GHZ proposal, three apertures are used 
to select the subset of particles emitted in the direction 
of the three analyzers. In our case all fields are emitted 
into the same spatial mode, so that one has to rely on 
beamsplitters to produce three outgoing beams. In such 
a scheme one cannot guarantee that photons are actu- 
ally registered by separate analyzers, as required by the 
GHZ argument. This was taken into account in the sim- 
ulations by discarding trajectories with a single detector 
being "hit" more than once. This way we could also dis- 
criminate against trajectories, in which more than three 
photons were detected (due to occasional diabatic tran- 
sitions within the dark space manifold). 

A comparison between the triple correlation function 
(L\((j)) £ 2 (<^>) Ls((f>)) , obtained by evaluating the quantum 
jump record of Monte Carlo simulations (with a decay 
rate k = 0.1T), and the result calculated for an ideal 
GHZ state [Eq. (21)], which would have been generated 
in the absence of cavity decay, is shown in Fig. 17. Equal 
angles 4> were chosen for all three detectors. Within the 
accuracy of the Monte Carlo simulation results, the quan- 
tum jump correlations precisely follow the cosine function 
predicted analytically for a stationary GHZ state. In par- 
ticular, perfect correlations are found for Y^k 4>k = Q and 
J2 k (f>k = 7r. The Monte Carlo wave-function simulations 
were repeated for different detector orientations. When- 
ever the sum over the angles was an integer multiple of 
7r, every trajectory with exactly three photons detected 
showed the correlations expected from a GHZ state, even 
when the field started leaking from the cavity well before 
the last photon was deposited to the cavity modes by the 
atom. 

It is remarkable that GHZ-type correlations are present 
between the decay photons even when the average occu- 
pation of the cavity mode is far below one. The observed 
correlations must be attributed to the dynamics of the 
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FIG. 17. Average triple correlations of photons decaying 
from the cavity, obtained from the quantum jump record of 
a Monte Carlo wave-function simulation. The adiabatic pas- 
sage parameters are chosen as in Fig. 14, but with a finite 
cavity decay rate k = 0.1F. All three analyzers are set at 
the same angle (0i = 02 = <f>z ), but equivalent results are 
obtained for different orientations. The solid line represents 
the correlations expected for an ideal GHZ state of the cavity 
field in the absence of decay. 

adiabatic passage and cannot be the property of a cer- 
tain quantum state of the cavity field [such as the case 
in Eq. (15)]. Preservation of coherence between subse- 
quent photon decays is clearly required and is provided 
by the Zeeman ground states, which are not affected by 
decay of the cavity field (see the time evolution of the 
atomic level population in Fig. 15 in Sec. V). The fact 
that entanglement of the "flying" photons is generated 
by the system dynamics distinguishes our method from 
most other proposals for the generation of GHZ corre- 
lated particles. Another notable feature of our scheme is 
that cavity dissipation does not lead to a loss of coher- 
ence, but rather is utilized for the birth of these flying 
"fields" with GHZ correlations into the external environ- 
ment. 

An important question is the dependence of the triple 
correlations on the cavity decay time. Simulating the 
adiabtic passage for different k, we found that the ap- 
pearance of perfect correlations persisted even if k was in- 
creased to a value of 10r, as long as only trajectories with 
exactly three decay photons were considered. However, 
the number of trajectories with five, seven, or more pho- 
tons generated by diabatic transitions during the atom's 
passage was found to grow with k. This is illustrated 
in Fig. 18, showing the distribution of the number of es- 
caped photons for four different cavity decay rates. Ob- 
viously, with more than three photons generated in the 
cavity, no triple correlations can be expected. While it is 
no problem to discard these trajectories in the numerical 
simulations, in an actual experiment these processes may 




FIG. 18. Distribution of the number of photons escaping 
during a single atomic transit for different cavity decay rates. 
Otherwise, the parameters are the same as in Fig. 14. For 
faster cavity decay, the probability for higher photon counts 
grows due to diabatic transitions within the dark state man- 
ifold. 



lead to wrong results due to finite detection efficiency. If, 
for example, two out of five emitted photons are not de- 
tected, the remaining photons are mistakenly registered 
as a three-photon event but would show no definite corre- 
lations. The probability for these events must therefore 
be minimized, which can be achieved by using a suit- 
able detuning and strength of the atom-cavity coupling 
to suppress diabatic transitions during the passage of the 
atom (as described in Sees. Ill and IV) and by chosing a 
small cavity decay rate. Since excited atomic states are 
occupied with an acceptably small probability only, the 
system is not compromised by atomic spontaneous emis- 
sion, which otherwise would also result in the generation 
of excess photons. 

In an actual experiment, one may take advantage 
of some simplifictions. For the experimental verifica- 
tion of the GHZ argument [cf. Eq. (22)], two of the 
three analyzers should be set at the same variable an- 
gle (</>i = 4>2 = 4>), while the other is always set at an 
angle fa — 0. Therefore, it is sufficient to use a pair of 
analyzers, with a half-wave plate only in front of one of 
them. Instead of using beamsplitters to route the decay 
photons, a cavity with two semitransparent output mir- 
rors may be used with an analyzer placed on cither side. 
The proposed setup is shown in Fig. 19. To avoid prob- 
lems due to detector dead time, events should only be 
considered if all the photons from the passage of a single 
atom arc registered by different photon counters. If the 
probability for the generation of more than three photons 
is negligible (which is the case, for example, for a detun- 
ing S s=s 0.6r and a small cavity decay rate k w O.ir, 
cf. Figs. 12 and 18), this last condition may be fulfilled 
by discarding runs with less than three distinct counters 
clicking. 
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FIG. 19. Simplified setup for a realistic experiment to 
probe GHZ correlations of photons emitted from a cavity with 
two output ports. The pair of detectors A analyzes the photon 
polarization in a rotated linear basis (given by the orientation 
of the wave plate), while detectors B use a fixed linear basis. 
The GHZ argument may be tested by observing the result 



of analyzer B if both detectors 
photon. 



and A y have registered a 



In this scheme, GHZ correlations could be detected 
by examining the result of the B counters when both A 
counters detect a photon. According to Eq. (21), the ex- 
pectation value for the outcome of the B measurement 
is then — cos(20), as the product of the A results is al- 
ways negative. Perfect correlations result for <f> = ir/2, 
when only the 5^-counter should receive the third pho- 
ton and for = 0, when the third photon should always 
be registered by B y . To maximize the probability for two 
photons being emitted through the mirror on the v4-side, 
and one through the B mirror, the A mirror transmission 
T should be twice that of the 5-mirror. In this case, the 
measurement record would show three separate photo- 
electric events in the required combination of analyzers 
in 44.4% of the atomic transits. 



VII. ATOM-PHOTON GHZ CORRELATIONS 

The method proposed in Sec. V for the synthesis of 
polarization-entangled states of the cavity field requires 
the initial generation of coherent superpositions of Zee- 
man ground-state sublevels. While certainly feasible, this 
requirement adds the need for additional microwave or 
laser pulse manipulation of the atom before the adiabatic 
passage. The state synthesis procedure would be greatly 
simplified if it were possible to generate entangled states 
starting from a single Zeeman sublevel, which is easily 
prepared by optical pumping. In this final section we will 
show that, by using a different Zeeman scheme, it is in- 
deed possible to generate a GHZ- type state with an atom 
initially pumped to the m g = ground state. However, 
in this case one of the entangled particles is the atom 
itself, which must be detected in order to observe GHZ 
correlations. 




FIG. 20. Level scheme of the (J 9 = 2) — » (J e = 1) transi- 
tion, showing the transitions relevant for an atom prepared in 
the m g — Zeeman sublevel (solid circle). The wide arrows 
indicate the 7r-polarized pump transitions. The predominant 
cavity-induced transitions are shown as solid arrows, the un- 
desired transitions are dotted. After the adiabatic passage, 
the atom is found in a superposition of the m g — —2 and the 
in g — +2-level (open circles), while the cavity contains two 
photons in a superposition of a + and a_ polarization. 



The atomic detection should be efficient, have a good 
time resolution, and a low background. These require- 
ments are met by using metastable atoms, which can 
be readily detected with secondary electron multipliers 
[48,49]. A suitable transition is ls 5 ( J = 2) — > 2p w (J = 
1) in argon, starting from the metastable state IS5, with 
a transition wavelength of 912.30 nm and a decay rate 
r/27r = 3.0 MHz. Figure 20 shows the corresponding 
J = 2 — > J = 1 level scheme. The adiabatic passage pro- 
cedure again uses a 7r-polarized pump field. The atom in 
this case must be prepared in the m = level, e.g., by 
pumping on the IS5 (J = 2) — ► 2p 8 (J = 2) transition, 
and the cavity modes are assumed to be empty. The low- 
est state of the dark space manifold that connects to this 
initial state is 



\Eo) 



-VT2g 2 (t)\g ,0,0) 
+20(t)fi(t)(|0_i,l,O) + | 5 i,O,l>) 
-fi 2 (t)(| fl _ 2 , 2,0) + |. 92 , 0,2)). 



(26) 



Due to the symmetry of the interaction Hamiltonian of 
Eq. (2) with respect to cr + /(T_ and positive/negative 
magnetic quantum numbers to, the adiabatic passage 
proceeds from the initial state to superposition states 
of positive Zeeman levels with cr_ photons and negative 
Zeeman levels with <r+ photons. In contrast to the case 
considered in Sec. V, the adiabatic passage leaves the 
atom in a state entangled with the cavity field: 



|fl0>®|0,0) =► 1.9-2,2,0) + \g 2 , 0,2). 



(27) 



As in Sec. IV, diabatic transitions may be sup- 
pressed and the system dynamics restricted to adia- 
batic changes of state (26) by using a finite cavity de- 
tuning 5. This is demonstrated by the results of a 
Monte Carlo wave-function simulation of the system us- 
ing the parameters for the above specified transition in 
argon. In Fig. 21, the occupation that is initially in level 



18 



CO 
_Q 
O 




laser 



SEM 



10 20 30 

Time t (in units of r 1 ) 

FIG. 21. Monte Carlo wave-function simulation of two-po- 
larization adiabatic passage in the case of a J = 2 - * J — 1 
transition. Atomic spontaneous emission is included, but cav- 
ity decay is neglected. The system evolution is given in terms 
of the occupation probability of the relevant basis states. Note 
that due to symmetry, the curves for the levels \g m , 0, m) and 
\g- m ,m,0) fall exactly on top of each other. The amplitudes 
of the Gaussian pulses for metastable argon were chosen as 
go = 30r (at t = 17r _1 ) and O = 50r (at t = 23r~ 1 ), with 
FWHM 10T -1 . The cavity detuning is S = 0.6T. 



\go,0,0) is transferred with equal probability to the lev- 
els \gi, 1, 0), \g~i, 0, 1) to end finally in a superposition 
of the states \g2, 2, 0), \g~ 2 , 0, 2), where for these calcula- 
tions cavity decay has been neglected. 

The final state in Eq. (27) should also display the cor- 
relations discussed by Grccnbcrgcr, Home, and Zeilingcr, 
which arc independent of the type of particle considered. 
The cavity decay photons are detected by using analyz- 
ers in front of the two cavity output ports as depicted in 
Fig. 19. The detection of the state-information carried by 
the atom requires a different strategy. By analogy with 
the optical case, it should involve the interference of pos- 
itive and negative Zeeman states and allow for a variable 
phase shift between both components. A possible exper- 
imental realization of such a detector is a Stern-Gcrlach 
analyzer with the magnetic field oriented perpendicular 
to the cavity axis (z axis), which is the spatial quanti- 
zation axis. The required phase shift is then obtained 
by rotating the magnetic field axis in the x-y plane by 
an angle (9/4. A sketch of the proposed setup is shown 
in Fig. 22. The energy eigenstates in the analyzer are 
obtained from the following operator involving angular 
momenta J± acting on the ground state 



J{6) = ^(e- ie/4 J + + e te / 4 J_), 



(28) 



with eigenvalues mj(9) = 0, ±1,±2. To establish GHZ 
correlations between the atom and the photons emitted 
from the cavity, one considers the operator with an eigen- 
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FIG. 22. Experimental setup for the detection of GHZ cor- 
relations between two cavity decay photons and the atomic 
state. The photons are detected as illustrated in Fig. 19. Af- 
ter the adiabatic passage, the atoms traverse a Stern-Gerlach 
detector with a magnetic field rotated by an angle (9/4. A 
grating then selects either the odd or the even orders of the 
Zeeman-split atomic states, which are then detected by a sec- 
ondary electron multiplier (SEM). 



value of 1 if the outcome of the Stern-Gerlach measure- 
ment is even [mj(9) = 0, ±2] and an eigenvalue of — 1 
if it is odd [mj(9) = ±1]. This operator may be ex- 
pressed as the difference of the projectors onto even and 
odd eigenstates of J (9): 



M(9) = P cvcn (9) - PoM = 
= e* \g. 2 )(g 2 \ + e^ 2 



+ e 



-iO/2 



9i){9- 



i>(ffi| + \go}(go\ 
ie \g 2 ){g-2\- (29) 



Note that the structure of M(9) is entirely equivalent to 
that of L{<j>) defined in Eq. (18), with the polarization- 
flip operators L± replaced by operators \g- m ) (gm\ chang- 
ing the sign of a Zeeman substate. From Eq. (19), 
the expectation value for triple correlations of the two 
cavity photons (using the results of Sec. VI) and the 
atom may be calculated, assuming the system is in state 
|GHZ')-(|g- 2 ,2,0} + |. 92 ,0,2))/V2, 

(GHZ'|L 1 (0 1 )i 2 (0 2 )M(0)|GHZ') = cos (fa + <j> 2 + 0) . 

(30) 

Thus, apart from the different origin of the third angle, 
Eq. (30) is exactly analogous to Eq. (21). Therefore, the 
conclusions of Sec. VI apply, in particular the discussion 
of perfect correlations, which occur when fa + fa + is 
an integer multiple of ir. 

Results for the actual experimental situation, with 
photons decaying from the cavity field before the adi- 
abatic passage is complete, must again be obtained by a 
Monte Carlo wave-function simulation of the full master 
equation. In addition to detecting the quantum jumps 
of the cavity field, the atomic ground state is probed at 
the end of each trajectory. An outcome for the mea- 
surement is chosen randomly, weighted with probabili- 
ties obtained from the final wave function, and the result 
is added to the classical record. Averaged triple corre- 
lations obtained as a function of the orientation of the 
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FIG. 23. Monte Carlo wave-function calculation of aver- 
aged triple correlations between two cavity decay photons and 
the atomic state after the adiabatic passage. The parameters 
are as in Fig. 21 and the cavity decay rate is k = 0.1F. The 
two photon analyzers are set at an angle <f>\ — <f>2 = and 
only the angle of the magnetic field 9 in the Stern-Gerlach 
detector is varied. The solid line represents the correlations 
obtained in the stationary state (|g-2, 2, 0) + \g2, 0, 2))/V2. 

magnetic field vector in the Stern-Gerlach detector are 
shown in Fig. 23. The cosine function of Eq. (30) is 
well reproduced, showing that the dynamic generation of 
GHZ correlations is also effective if the atomic degrees of 
freedom are entangled with the photons. 

From an experimental point of view, GHZ correlations 
between an atom and two photons may be easier to de- 
tect than their three-photon counterpart. With only two 
photons escaping from the cavity, the fraction of usable 
detector responses would increase to 50%. In addition, 
detection of the atom provides a trigger signal that can 
be used to discriminate against background noise in the 
photon counters. Clearly another advantage is that there 
is no need for the preparation of a superposition of Zee- 
man ground states prior to the adiabatic passage. 

VIII. CONCLUSIONS 

We have considered a model of adiabatic passage in a 
strongly coupled atom-cavity system exposed to a clas- 
sical pump field. With such a setup one may synthesize 
nonclassical photon states by mapping Zeeman ground- 
state coherences onto the cavity mode. The original pro- 
posal presented in Refs. [20,23] included only a single 
cavity mode. We have extended that scheme to include 
two degenerate cavity modes with orthogonal polariza- 
tion, the case relevant to Fabry-Pcrot resonators used 
in the laboratory. By analytical calculations as well as 
quantum Monte Carlo simulations we have shown that 
atomic passage through the system does not lead to pure 



number states of the two-mode cavity field, as long as 
atom, cavity, and external pump field are resonant with 
each other. We could identify diabatic transitions to de- 
generate dressed levels with higher photon number as the 
source of this effect. 

A central result of the paper is that the adiabatic pas- 
sage scheme can be modified to work even in the two- 
polarization case by introducing a finite cavity detuning. 
This way diabatic transitions are suppressed and states 
with a definite photon number of the cavity field may be 
efficiently generated. Several applications of the method 
were discussed. As a first test of our scheme, we cal- 
culated the synthesis of a three photon Fock state of a 
single polarization of the cavity field and obtained an effi- 
ciency above 99% based on parameters actually achieved 
in cavity-QED experiments. 

The system under study offers intriguing possibilities 
beyond the generation of number states and their super- 
position for a single mode. An important property of 
the atomic transition scheme we employ is its symmetry 
under a change of sign of (circular) polarization and mag- 
netic quantum number. Taking advantage of this feature 
allows one to create polarization entangled states of the 
cavity field. As a particularly relevant example we pre- 
sented a proposal for the synthesis of a so-called GHZ 
state consisting of three photons maximally entangled in 
polarization. It can be used to test EPR correlations 
without having to resort to Bell-type inequalities. 

To probe correlations between the photons generated 
by our scheme, we make use of the inevitable presence of 
cavity dissipation. The decay photons emitted from the 
cavity are sent to linear polarization analyzers of variable 
orientation. For certain combinations of detector orien- 
tations, perfect correlations between the photons may 
be observed. Most remarkably these correlations per- 
sist even if cavity decay and state synthesis occur on the 
same time-scale. 

One should note a difference between the states syn- 
thesized in the cavity and the states discussed by GHZ 
[3]. In the latter case, each particle occupies a distinct 
mode, while in our cavity the photons are, in princi- 
ple, not distinguishable. Nonetheless, as demonstrated 
in Figs. 17 and 23, the observed correlations for a par- 
ticular subensemble would be precisely as for the "ideal" 
GHZ state. Qualitatively, this correspondence is due to 
the temporal evolution of the adiabatic passage and the 
stochastic nature of cavity decay with photons usually 
emerging from the resonator "one by one." This effect 
could be enhanced through manipulation of the adiabatic 
passage process: combining suitable shifts of the atomic 
states with multipeaked interaction profiles g(t) and Cl(t), 
the adiabatic passage would proceed in well-separated 
stages with exactly one photon emitted per stage. With 
the photons arriving at the detectors at distinct (prede- 
termined) times, the observed correlations generated by 
the adiabatic passage should correspond exactly to those 
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predicted by GHZ without the need for post selection of 
a particular subensemble. 

A practical concern for the current scheme is the lim- 
ited time resolution of the detectors, which may not be 
sufficient to separate the decay events. In this case one 
might simply register only events in which each of the 
three photons is actually detected by a separate detector. 
Note that, in the three particle interferometer orginally 
proposed by GHZ [4] , the correct direction and energy of 
the particles must also be enforced by apertures around 
the source and filters in front of the detectors. 

We did not discuss the influence of finite detection effi- 
ciency on the correlations observed. However, as long as 
the number of photons generated does not exceed three 
(a condition that, as we have demonstrated, may be ful- 
filled by chosing a sufficiently large cavity detuning or co- 
herent coupling go), missing a photon decaying from the 
cavity would only lower the coincidence rate observed, 
but would not lead to modified three photon correlations. 
The reason is that events with less than three detected 
photons can easily be discarded in a coincidence mea- 
surement. This argument is sometimes called the fair 
sampling assumption [40]. 

The proposed method for the synthesis of entangled 
states is not restricted to the specific cases discussed in 
the paper. More complicated entangled states of the cav- 
ity field may be generated as well. Examples include su- 
perposition states with different photon number in the 
two polarization modes (e.g., |3,0) + |0, 2)) or entangled 
states with polarization distributed in a more general way 
(e.g., |4, 2) + |1, 3)), as would be useful in quantum infor- 
mation networks with error correction [26]. 

Maximally entangled states with more than three par- 
ticles, as discussed by Mermin [50], are of particular in- 
terest, since they provide an exponentially stronger test 
of local realism. These states may be generated with 
our scheme without any modification, if an atomic tran- 
sition with the corresponding angular momentum quan- 
tum number is used. In cesium, a four particle gener- 
alization of the GHZ state may be synthesized with the 
F g = 4 — > F e = 4 hyperfine component of the D2 line. 
However, with increasing number of particles, detection 
of the correlations between them becomes considerably 
more challenging, since the coincidence rate would be 
severely reduced. 

The possibilities for the generation of entanglement 
in our system are further expanded by including atom- 
photon correlations in the considerations. A correspond- 
ing experiment would offer at least two additional bene- 
fits: instead of a coherent superposition of Zeeman sub- 
levels, only a single atomic substate is needed as an initial 
condition for the adiabatic passage, so that the prepara- 
tion of the atomic state would require only standard opti- 
cal pumping. Furthermore, detection of the atom would 
provide a gate for the photon detectors, increasing the 
signal to noise ratio. 



In this paper we have restricted our attention to the 
case of a dilute atomic beam passing through a cavity. 
In view of experimental progress in the field [32,33], one 
may also consider atoms moving through the cavity at 
a low enough speed to be considered stationary on the 
timescale set by the coherent interaction, or even trapped 
atoms or ions in a cavity. The method presented here 
could also be applied to the case of stationary atoms in a 
cavity, if the dynamics that in ordinary adiabatic passage 
is induced by the atom's motion, is generated instead by 
suitably timed pulses. The interaction with the two cav- 
ity modes may be controlled by coupling the atomic lev- 
els with a Raman transition using the cavity field and an 
auxiliary classical field. As suggested in Refs. [22,25,26], 
the Raman coupling strength may then be adjusted by 
changing the intensity of the auxiliary field. 
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